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_N  A  V  *v- 


SOME  RELIABILITY  APPLICATIONS  OF  THE  VARIABILITY  ORDERING 

by 

Sheldon  M.  Ross  and  Zvi  Schechner 

1.  A  VARIABILITY  RESULT 

If  X^  and  X2  are  random  variables  having  respective  distributions 

F^  and  ?2  »  then  we  say  that  X^ 

Xj)  or  equivalently  that  <_ 

v 

00 

J  f(x)dF1<x)  < 

—00 

for  all  increasing  convex  functions  f  .  Some  easily  derived  properties  of 
this  ordering  are 

1.  F.  <  F2  if  and  only  if 

A  v 

00 

j  f1(x)dx 

a  a 

where  f ^  ■  1  -  F^  . 

2.  If  F^  <_  G^  ,  i  -  1,2  ,  then  F^  *  Fj  <,  G^  *  ®2  w*1*re  *  denotes 

v  v 

convolution. 

He  will  now  present  a  theorem  concerning  this  ordering  and  in  Sections 
2  and  3  apply  it  to  branching  processes  and  shock  models.  Other  applications 
of  the  variability  ordering  to  population  survival  models  (Section  4)  and  to 
Poisson  shock  models  (Section  5)  will,  then  be  presented. 


OB 

J  f(x)dF2(x) 


<  X«  (read  X.  is  less  variable  than 
v  *  1 


Thgorgi  1: 

Let  X^^,  . . .  be  a  sequence  of  nonnegative  Independent  and  identically 

distributed  random  variables  and  similarly  .  Let  N  and  M  be 

integer  valued  nonnegative  random  variables  that  are  independent  of  the  {X^} 
and  {Yj}  sequences.  Then 


N  M 

l  >  l  • 

V  V  1-1  v  1-1 


Proof: 


We  will  first  show  that 


l 

1-1 


M 

Xi-  K 
1  v  1-1 


Let  h  denote  an  increasing  convex  function.  To  prove  the  above  we  must 
show  that 


<»  *£8  *i)]  i  {4  xJ]  • 

Since  H  >  M  ,  and  they  are  Independent  of  the  X^  ,  the  above  will  follow 
v 

If  we  can  show  that  the  function  g(n)  ,  defined  by 

g(n)  -  E[h(X^  +  ...  +  Xn)] 


is  an  increasing  convex  function  of  n  .  As  It  is  clearly  increasing  since 
h  is  and  each  X^  is  nonnegative  It  remains  to  show  that  g  is  convex, 
or,  equivalently,  that 


(2) 


g(n  +  1)  -  g(n)  is  Increasing  in  n 


To  prove  this  let  sn  "  1  ^  »  aiM*  note  that 


g(n  +  1)  -  g(n)  -  E[h(Sn  +  X^)  -  h(Sn>] 


E[h(Sn  +  X^)  -  h(S„)  |  Sn  -  t]  -  E[h(t  +  X^,)  -  h(t)] 


n'  ■  n 


£(t)  (say). 


As  h  la  convex,  it  follows  that  f (t)  is  increasing  in  t  .  Also, 

S  increases  in  n  ,  we  see  that  E{f(S  )]  Increases  in  n  .  But 

n  n 


E[f(Sn)]  -  g(n  +  1)  -  g(n) 


and  thus  (2)  and  (1)  are  satisfied. 
We  have  thus  proven  that 


?  r 

I  \  >  Z  *t 

1  1  v  1  1 


and  the  proof  will  be  completed  by  showing  that 


M  M 

Z  xi  >  Z  *t 

l  Si  1 


or,  equivalently,  that  for  Increasing,  convex  h 


But 


*  h(Z  *i)  h)  ■ 


by  Independence 


■  m 

since  I  X,  >.  X  Y. 
1  1  v  1  1 


and  the  result  follows  by  taking  expectations  of  both  sides  of  the  above 


!♦  A  BRANCHING  PROCESS  APPLICATION 

Consider  two  Galton  Watson  branching  processes  in  which  individuals 
at  the  end  of  their  lifetime  give  birth  to  a  random  number  of  offspring. 
Let  »  j  il  »  n  >  0  denote  the  number  of  offspring  of  the  jth 

individual  of  the  nC^  generation  in  the  i^1  branching  process ,  i  ■  1  ,  2 
Suppose  that  the  random  variables  X^  i  3  >1  »  n  >  0  are  independent 
for  i  -  1  ,  2  and  have  a  distribution  not  depending  on  j  .  In  addition 
suppose  that 

xjj*  L  £or  all  n  ,  j  . 

Let  ,1*1,2  denote  the  size  of  the  n^  generation  of  the  l**1 

process. 

Proposition  2: 

If  zjS1*  -  1  ,  i  -  1  ,  2  ,  then  Z*1*  >  zf2)  for  all  n  . 
u  n  “  u 


Proof: 

The  proof  is  by  induction  on  n  .  As  it  is  true  for  n  *  0  , 
assume  it  for  n  .  Now, 


ZU> 

Z(1)  -  j  - 

n+1  fa  J»n 


and  so  the  result  follows  from  Theorem  1. | j 


Let  P{X  -  1}  -  a  ,  P{X  -  0}  -  (1  -  a)  -  ,  P{X  -  2}  -  M-y° 

and  let  Y  be  e  nonnegative.  Integer  valued,  and  such  that  P{Y  -  1}  <_  a 

and  E[Y]  -  M  .  If  a<M<2-ct,  then  X  <  Y  . 

v 

Proof: 

We  mist  show  that 


l  P{X  >_  1}  <  l  P{Y  >.  1}  ,  n  -  1,2,  ...  . 
i-n+1  i-n+1 


As  E[X]  -  E[Y]  »  M  .  this  Is  equivalent  to 


l  P<X  i  1}  ^  l  P{Y  >  1}  ,  n  »  1,2,  .... 

1-1  1-1 

When  n  -  1  ,  the  above  reduces  to  P{X  -  0}  <_  P{Y  -  0)  .  This  follows 
since,  as  P{Y  -  1)  <  P{X  -  1}  ,  if  P{Y  -  0}  <  P{X  -  0}  then  it  would 
not  be  possible  for  E[Y]  to  equal  E[X]  .  When  n  >  1  ,  the  above  is 
equivalent  to 

n 

M  >  l  P{Y  >  1} 

1-1 


which  follows  since  E[Y]  -  M  .|j 


A  SHOCK  MODEL  APPLICATION 


Suppose  that  shocks  occur  in  accordance  with  a  renewal  process  having 
interarrival  distribution  G  and  mean  .  Each  shock  gives  rise  to  a 
nonnegative  random  damage  which,  independent  of  all  else,  has  probability 
distribution  F  .  The  damages  are  assumed  to  be  additive  and  we  let 
D(t)  denote  the  damage  at  time  t  .  That  is, 

M(t> 

D(t)  -  I  X 
i-1  1 

where  X^  is  the  damage  of  the  i^1  shock  and  M(t)  is  the  number  of 
shocks  by  t  .  The  system  is  assumed  to  fail  the  first  time  that  D(t) 
exceeds  some  constant  c  .  That  is,  the  system  falls  at  time  Tf  G  where 

Tj,  G  -  min  {t  :  D(t)  >  c)  . 

We  will  obtain  a  variability  result  about  T?  G  when  both  F  and  G 
are  NBUE  distributions,  where  a  distribution  of  a  nonnegative  random 
variable  X  is  said  to  be  NBUE  (new  better  than  used  in  expectation)  if 

E[X  -  t  |  X>  t]  <  E[X]  for  all  t  >.  0  . 


Letting 

N(c)  ■  max  {n  :  X^  +  ...  +  X^  <_  c}  . 

Then  the  system  will  fail  at  the  time  of  the  N(c)  +  1  shock. 

Lemma  5: 

If  F  is  NBUE,  then 

N(c)  +  1  <  N*(c)  +  1 
v 


•-.iv 


aaHEaaaSL-atiifti^sj 


MM2 


»_**%.-•%  »*  *.  ’-.  \  -M-." -. 


where  H  (c)  le  e  Poisson  random  variable  with  mean  c/yp  where  u^,  ■  E[X] 


Proof : 

As  N(c)  is  just  the  number  of  renewals  by  time  c  of  a  renewal  process 
Whose  interarrival  distribution  is  NBUE  with  mean  yp  ,  the  result  follows 
from  Theorem  3.17  on  page  173  of  [1].| | 


Proposition  6: 


If  F  and  G  ere  both  NBUE  distributions,  then 


T  <  T 

tf.o  -  \.E2 


where  E^  and  ere  exponential  random  variables  having  the  same  means 
as  F  and  G  respectively. 


Proof: 


We  can  express  Tp  G  by 


N(c)+1 


r,G  "  **■ 

i"l 


where  the  Y^  ,  i  >_  1  ,  are  the  Interarrival  times  between  successive  shocks. 
They  ere  thus  Independent  and  have  distribution  G  .  Now.lt  is  well  known 
that  an  NBUE  distribution  G  is  less  variable  than  an  exponential  distribu¬ 
tion  with  the  same  mean  and  so 


Y.  <  e .  when  e.  is  exponential  with  mean  y_  . 
l—ii  G 


The  result  now  follows  from 


3  end  Theorem  1. 


As 


E1’E2 


N  <c>+1  * 

I  e.  ,  it  follows  upon  conditioning  on  N  (c)  that 
i-1  1 


where  G  (x)  is  the  goal  distribution  with  paraoeters  n  and  l/y_ 

(its  mean  is  ny^)  .  Also,  if  F  and  G  are  NBUE,  then  from  Proposition  6 

all  of  the  somts  of  T  _  are  no  greater  than  the  corresponding  moaents 

r  »b 

of  I  .  ,  For  instance, 

*1»*2 

-  E[(N*(c)  +  l)y  ]  -  (c/y_  +  l)y_  . 


e[i.  „i  « ,  1 


Consider  a  population  of  m  individuals  each  of  whom  is  required  to  spend 


exactly  one  time  unit  out  In  the  field.  For  each  day  i  ,  i  -  1,2,  ...,  m  , 
there  is  a  random  variable  Y^  which  represents  the  probability  that  an  indi¬ 
vidual  out  on  the  field  on  day  i  will  survive.  That  is,  given  Y^  ,  each 
individual  sent  out  on  day  i  will  independently  survive  with  probability 
Y^  .  The  Y^  ,  i  ■  1,  . ..,  m  are  assumed  to  be  independent  and  Identically 
distributed  random  variables  for  which  P{0  <_  Y^  <_ 1}  ■  1  . 

A  strategy  for  the  population  is  a  positive  integer  valued  vector  n  - 

k 

(n^,  . . . ,  n^)  ,  k  <_  m  ,  £  n^  ■  m  ,  with  the  interpretation  that  n^  individuals 
are  sent  out  on  day  i  ,  i  ■  1,  . ...»  k  .  Let  N(n)  denote  the  number  of  indi¬ 
viduals  that  survive  under  strategy  n  . 

Proposition  7: 

H(n)  >1(1,1,  ...,  1)  . 
v 


Proof: 

As  the  variability  ordering  is  closed  under  convolution,  it  clearly 
suffices  to  prove  that  if  n  individuals  are  sent  out  on  a  given  day,  then 
the  nuaber  of  survivors  is  more  variable  than  it  would  be  if  the  n  indi¬ 
viduals  were  sent  out  on  separate  days.  Hence,  we  must  show  that  for  any 
convex  function  f 

*[£  f<i)(J)Y4<l  -  I)*^J  >  f(i)(J)[H(Y)]i[l  -  E(Y) . 

However,  it  follows  from  the  following  lemma  that  the  function 


g(p)  =  I  f(i)(V(i  -  p)1 

i-0  'x/ 


is  convex  in  p  and  ao  the  above  follows  from  Jensen's  Inequality.! 


Let  *  *  *  be  independent  Bernoulli  randoa  variables  with 

P{X^  -  1}  ■  p  .  If  f  is  a  convex  function,  then. for  any  constant  c 


‘[‘(I +  c) 


is  a  convex  function  of  p  . 


Proof: 


The  proof  is  by  induction.  When  n  ■  1  we  wust  show  that  pf(c  +  1)  + 


(1  -  p)f(c)  is  convex  in  p  ,  which  is  1— diets.  Assusing  the  result  for 

n  -  1  we  then  have,  upon  conditioning  on  X  , 

n 

*[*(! x* + c)]  ■  4(T  xi + c + x)] + <i  ■  4nl  x* + c)]  ■ 


Hence  if  we  let 


gc<P) 


'  '[<?  x‘ +  c)] 


ws  suet  show  that 


h(p)  =  Pg^jCp)  +  Cl  -  p)gc(p) 


is  convex.  Differentiation  with  respect  to  p  yields 


h"<p)  -  2(g^.1(p)  -  g;<p)>  ♦  pg^Cp)  +  a  -  p)g"<p> 


J-'?- 


13 


Mow,  gg'Cp)  and  g^J^Cp)  are  nonnegative  by  the  Induction  hypothesis 
Also,  aa  £'  Is  Increasing  In  x  ,  we  see  that 

*:<«->  -  ^(T  +  *)] 

is  Increasing  in  a  and  thus  g'  .  (p)  -  g'  (p)  >0  .  Hence, 

CTl  c  “ 

h"(p)  >0 

which  proves  the  result.  1 1 
Renark: 

Since  E[N(n) ]  -  mE[Y]  for  all  strategies  n  ,  it  follows,  as  in 
Corollary  3,  that  P{N(g)  ■  0}  Is  minimised  by  the  strategy  (1,1,  .. 


14 


5.  A  POISSOH  SHOCK  MODEL 


Suppose  thet  unobserved  shocks  hit  s  device  In  accordance  with  a  Poisson 
process  having  rate  X  .  At  some  tine  the  device  is  inspected  and  If  n 
shocks  have  occurred  by  that  time  a  cost  f(n)  Is  Incurred,  where  f  is  an 
Increasing,  convex  function.  There  are  two  inspection  plans  that  can  be 
employed — one  of  which  inspects  the  device  after  a  random  time  X  and  the 
other  after  a  random  time  Y  ,  and  we  are  interested  in  determining  which 
plan  leads  to  a  smaller  expected  cost. 

Proposition  9  shows  that  if  X  is  more  variable  than  Y  ,  then  it  leads 
to  a  greater  expected  coat. 


Proposition  9: 


Let  {N(t>  ,  t  >_  0}  be  a  Poisson  process  with  rate  X  and  let  X  and 
Y  be  nonnegative  random  variables  that  are  independent  of  the  Poisson 
process.  Then 


X  >  Y  e  N(X)  1  H(Y)  . 


Proof: 


Let  f  be  a  convex,  increasing  function  and  suppose,  without  loss  of 

/ 

generality,  that  f(o)  •  0  .  Let 


A(n)  -  f(n  +  1)  -  f(n) 


H(t)  -  l  f(n)(Xt)n/n! 
n«0 


and  note  that  A(n)  L  0  and  Increasing  in  n  by  assumption.  How, 


E[f(N(X))3 


«"XtH(t)dPx(t) 


where  ?x  is  the  distribution  function  of  X  .  Hence, it  suffices  to  show 
Chet  g(t)  ■  e_tH(t)  is  s  convex  increasing  function  of  t  .  How, 

g'(t)  -  -Ae"XtH(t)  +  e“XtH'(t) 

■  e“Xt  l  [f(n  +  1)  -  f  (n)]An+1tn/n! 
n«0 

-  AE[A(N(t))]  . 

Now,  A  >_  0  and  so  g'  >.  0  .  Also,  as  H(t)  is  stochastically  increasing  in 
t  and  A  is  increasing,  it  follows  that  E[A(N(t)) ]  is  increasing  in  t 
and  so  g  is  convex.  |  j 

1— rk: 


Proposition  is  not  true  for  a  general  renewal  process.  It  is  clearly 

not  true  fox  a  deterministic  renewal  process.  (If  all  Interarrivals  equal  1 

and  X  uniform  (1.9,2. 1),  Y  5  2  ,  then  X  >.  Y  but  N(X)  <  N(Y».  Even 

v  st 

if  the  lnterarrival  distribution  has  a  decreasing  failure  rate, Proposition  9 
need  not  be  true.  For  a  counterexample ,  suppose  the  lnterarrival  distribution 

la 

-A.x  -A.x 

F(x)  -  p(l  -  e  1  )  +  (1  -  p)(l  -  a  ^  )  ,  0  <  p  <  1  ,  Ax  ^  A2  . 

To  compute  E[N(t)]  ,  Imagine  that  at  each  renewal  a  coin,  having  probability 
p  of  landing  heads,  is  flipped.  If  head  appears,  the  next  lnterarrival  is 
exponential  with  rata  A^  and  if  tails,  it  is  exponential  with  rate  A2  • 
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I[M(Xe)]  -  A-  A  J  ^dx  +  B 
1-e 

-a[i.+  £  c.~  ■-.«»]♦ 

Hence,  E[N(X_)]  decreases  In  e  when  e  is  near  1 

6 

whenever  0  <  Cj  <  <  1  . 


B  . 


However,  X  >  X 
C1  v  62 
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